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Inflow and outflow from the accretion disc of the microquasar 
SS433: UKIRT spectroscopy 
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ABSTRACT 

A succession of near-IR spectroscopic observations, taken nightly throughout an entire 
cycle of SS 433 's orbit, reveal (i) the persistent signature of SS 433 's accretion disc, having a 
rotation speed of ~ 500 km (ii) the presence of the circumbinary disc recently discovered 
at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster 
outflow than has previously been measured for the disc wind. From these, we find a much 
faster accretion disc wind than has noted before, with a terminal velocity of ^1500 kms~^. 
The increased wind terminal velocity results in a mass-loss rate of ~10~^ Mq yr~^. These, 
together with the newly (upwardly) determined masses for the components of the SS 433 
system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington 
flows. Our observations imply that the size of the companion star is comparable with the 
semi-minor axis of the orbit which is given by Vl — 40 Rq , where e is the eccentricity. Our 
relatively high spectral resolution at these near-IR wavelengths has enabled us to deconstruct 
the different components that comprise the Brackett-7 line in this binary system, and their 
physical origins. With this line dominated throughout our series of observations by the disc 
wind, and the accretion disc itself being only a minority (~15 per cent) contribution, we 
caution against use of the unresolved Brackett-7 line intensity as an "accretion signature" in 
X-ray binaries or microquasars in any quantitative way. 



Key words: accretion, accretion discs 
binaries: spectroscopic 



stars: individual: SS 433 - stars: winds, outflows 



1 INTRODUCTION 

SS433, named in 1977 as Stephenson and Sanduleak's 433rd en- 
' try in their catalogue of Ka emitters, was first r ecognised as an 
X-ray and radio source bv lClark&Murdinl ( ll978h . possibly assso- 
ciated with the supernova remnant W50. It became famous as the 
first known relativistic jet source in the Galaxy, and is the only one 
kn own with bary onic content in its jets (see comprehensive review 
bv lFabrikall200i) . Its jets are believed to emerge from an accretion 
disc, but direct evidence for this has thus far been scant. 

The mass of t he system has been el usive over the last 30 
years. For example, iHillwig & GiesI (l2008h . from observations of 
absorption lines, reported masses of ~ 12 and ~ 4 M© for the 
donor star and the comp act object, respectively. On the other hand, 
iD'Odorico et al.l (Il99lh had previously suggested that the masses 
of the components were rather low, ^3.2 Mc?) and ^ 0.8 Mct) fo r 
the donor and the compact object, respectively. iLopez et al.l ( [200^ , 
looking at the eclipses seen in X-ray spectra, suggested masses of 
^ 35 and ^ 20 M© for the donor and the compact object, respec- 
tively. Recent analysis of the stationary Ha line has shown, from 
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dynamical considerations, that the enclosed mass for the whole sys- 
tem is quite large, around 40 Mq , and that the mass for the compact 
object plus accretion disc attains 16 Mq, therefore also implying 
that the id entity of the compac t object is a rather massive stellar 
black hole (iBlundell et aDl2008h . 

At near-infrared wavelengths, SS433 is characterise d by a 
red continuum and bright emission lines (a Thom pson et aD[T979l : 
iMcAlarv & McLarenlll98Ql) . As at optical and X-ray wavelengths, 
these emission lines can be divided into two groups: (a) stationary 
lines, although highly variable in strength and profile shape; (b) 
moving lines, thought to originate in two oppositely-directed rela- 
tivistic jets of baryonic content (moving with speed v ^ 0.26 c). 
The stationary emission lines are thought to be produced in the ac- 
cretion flow and in an expanding, geometrically thick environment 
fed fro m the stellar wind and the supe r-Eddington accretion disc 
outflow (iGies et al.|[2002l : lFabrikall2004l) . 

The system shows three main periodicitie s: the binary's orbital 
motio n, with a period of about 13.1 days ( Cr ampton & Hutchingsl 
Il98lh : and the jet axis's precession and nutation, with periods of 
about 162 and 6 days, respectively. A configuration where the jet 
axis undergoes a precession cycle in 162 days, referred as the kine- 
matical model, was proposed to describe the motion of the lines 
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(iFabian & Reesll979l;lMilgromll979h . In order to account for the 6- 
day nodding cycle jKatz et alj (Il982h proposed a dynamical model 
where the companion exerts a gravitational torque on the disc pe- 
riphery. 

Near-IR light can escape high opacity and dusty environments 
more easily than Ha photons can, making it possible to detect 
heavily obscured line-emitting regions. We carried out a line de- 
blending procedure based on fitting Gaussian profiles to our near- 
IR spectra, choosing Brackett-7 (hereafter Br7) at A 2.165 /xm as 
the most suitable emission line to model. 

The Br7 recombination line has been extensively used in the 
past to diagnose the presence of accretion discs an d outflows (e.g., 
IBandvopadhvav et al.l [19971 : IShahbaz et al]ll999h. In the case of 
cataclysmic variables (CVs), Dhillon & MarshI ( 1 995h argued that 
since the Brackett and He I emission lines are so strong in emission 
and also so broad that they must originate in the accretion disc. In 
low-mass X-ray binaries (LMXBs) the Br7 line has been detected 
showing a double-peaked pr ofile characteristic of accretion discs. 
IBandvopadhvav et aP (Il997h used the distance between the peaks 
to calculate properties of the orbits in LMXBs, assuming that the 
compact object accretes via Roche-lobe overflow. 

In this work, we present unprecedentedly high signal-to-noise 
mid-resolution near-IR spectra of SS 433 from UKIRT that has en- 
abled us to identify the accretion disc and its outflow. In Section |2l 
we begin by giving a brief description of the UKIRT observations 
and the data reduction. In Section[3]we discuss the results of decon- 
structing the Br7 emission line in order to identify and quantify the 
different components present in emission in the system throughout 
an orbital period. Finally we discuss and give concluding remarks 
on this work in Sections|4]and[5] 



2 OBSERVATIONS AND DATA REDUCTION 

We observed SS 433 every nigh t over an entire orbital cycl e with 
the UKIRT UIST spectrograph jRamsav Howat et all |2004|) from 
2006 August 17 to August 29, during which time the precession 
phase varied from 0.40 to 0.47 respectively. We use the conven- 
tion in which orbital phase {(j)ovh) zero is when the don or star is 
eclipsing out the compact object (iGoranskii et al.lll998l) . Preces- 
sional phases (^pre) are based on the ephemeris reported in lPabrikal 
(I2OO4I) . where precession phase zero is when the jet lines are max- 
imally separated hence the inclination of the jet axis with our line- 
of-sight attains a minimum, i.e., it corresponds to maximum expo- 
sure of the accretion disc to the observer. 

UIST is a 1-5 /xm imager- spectrometer with a 1024 x 1024 
InSb array. The plate scale in spectroscopic mode is 0^.^12 pixel" ^. 
We used a slit of width 0^.^24. The spectral resolution R — A/ A A 
is about 4000 for the J-band and 3900 for the iT-band. The top 
of the slit was pointing to the east and centred on SS 433 's optical 
position during acquisition. There were no problems in acquiring 
the target since SS 433 is a very bright red object with K mag. 
Our data consist of 13 nights of iT-band high-resolution spectra 
and one night, at primary minimum, covering the complete near-IR 
window (J, H and K -bands). 

The observing strategy consisted of taking about forty frames 
with exposure times of 45 s each, following a ABBA nodding pat- 
tern. The pixel-to-pixel response of the array was corrected using 
the normalised, reduced flat field. Bad pixels, correction of detec- 
tor defects and cosmic-ray removal were performed using standard 
median filtering techniques. Frequent observations of a standard 
star (BS6697) and the arc lamp were performed during each night 
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Figure 1. Individual Br7 stationary profiles plotted as a function of wave- 
length and modified Julian date of observation. Their continua are nor- 
malised to unity and aligned with the mean day of their observation. Each 
spectrum is labelled with its orbital phase on the right hand side. The dashed 
line corresponds to the systemic velocity measured from the Mg ii radial ve- 
locity curve (Perez et al. in prep). 



to provide an accurate estimate of the wavelength axis. Wavelength 
calibration was carried out by fitting a high-order polynomial to 
the Argon lamp arc spectra, taken each night. In order to remove 
the sky contribution, pairs of nodded frames were subtracted. Stan- 
dard star division was applied to each frame in order to flux cal- 
ibrate the data. The intrinsic photospheric features were carefully 
removed from the standard star by interpolating or fitting a Voigt 
profile when needed. Great care was taken to cancel the telluric 
features, especially the CO2 absorption near 2.05 /im. Finally, the 
spectra were extracted using Home's optimal extraction algorithm 
(lHornelll986l) and then combined to produce an average spectrum 
for each night with an exposure time of about 2 ks giving a signal- 
to-noise ~ 400 in the Br7 line. 

All data reduction and analysis were carried out using the Perl 
Data Language^ The reduced spectra centred on the heliocentric 
Br7 line are displayed in Fig. [T] from which it is immediately ap- 
parent that there is considerable structure in this emission line, and 
considerable variation of this structure. We deconstruct these struc- 
tures in the next section. 



3 DECONSTRUCTING THE BRACKETT-7 PROFILE 

Many other au thors have studied t he stationary Ha line profile 
in detail (e.g., iFalomo etaDll987l : iGies et all l2002h . albeit over 
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Figure 2. Upper panel: Example of 617 stationary emission lines observed 
at orbital phase (j)orh = 0.27; the top x-axis corresponds to heliocentric 
speed in units of kms~^. Lower panel: Tracks of the centroids of the 
Gaussian components fitted to each of our spectra. The modified Julian date 
(MJD) increases vertically and the tick mark heights are proportional to the 
FWHM of each component (see inset in the bottom left corner). 



rather shorter observing runs. Decomposition of a line profile as 
the sum of Gaussian components is a technique widely used to ex- 
tract information from different parcels of gas in a line emitting 
region. For example, in the case of Seyfert galaxies it allows one 
to di stinguish the na r row line region from the broad line region 
(e.g.. IHo et al.lll997h . iBhindell et all (120085 decomposed the sta- 
tionary Ka line (observed during a quiescent period in SS433's 
behaviour) into primarily 3 components: one broad component 
(FWHM ~ 700 km s~^) whose width was observed to decrease with 
precessional phase (i.e., as the jets become more in the plane of the 
sky) identified as the accretion disc wind, and two narrower, red- 



and blue- shifted components, but stationary in wavelength, being 
radiated from a glowing circumbinary ring. 

The stationary Br7 emission lines show a much more comple x 
profile than the quiescent Ha line studied bv lBlundell et 
After trying with different numbers of Gaussian components we 
came to the realisation that up to six components were needed to 
account for the complexity of the Br7 profile shape. Fig. [2 (upper 
panel) shows an example of a fitted Br7 profile; in this case six 
Gaussians were used in the fit. 

The presence of a P Cygni feature in the stationary lines 
has been noted by several authors at certain precessional phases 
dCrampton & Hutchingsl Il98ll : iFilippenko etaP Il988l : iGies et all 
l2002h . An absorption feature in the blue wing of the line pro- 
file would indeed complicate the analysis and it would have to 
be taken into account using models of outflowing winds (e.g., 
ICastor & LamersI 19791) . However, we have found no clear evidence 
of the presence of such absorption feature at the epochs at which 
we observed SS 433 (see e.g., top panel of Fig. [J]). 

Fig. [2] (lower panel) shows the components of the stationary 
Br7 line as a function of time. It is easy to see that the Br7 com- 
plex can be decomposed in three main constituents: a very broad 
wind component present at all times in our data- set and two sets 
of narrower pairs. The broad wind and both narrower pairs show 
a mean velocity close to the systemic speed Vb — 150 kms~^ 
(see lower panel in Fig.O. This measurement of the systemic ve- 
locity was obtained from the radial velocity curve of the Mg II A 
2.404 /im stationary line (Perez et al. in prep). Although this value 
seems to fit well with all the stationary emission lines in our near-IR 
spectra, it is very di fferent compared to other previo us observations 
(e.g., 27 kms" -^ bylCrampton & Hutchingsl (Il98ll), 73 kms ~^ by 
lHillwig&Giesl(l2008h and -44 kms'^bv lGTes et aLl ( l2002h ). 

The broad components plotted in Fig. |2] show FWHMs from 
1300 up to 1500 km s~^. The presence of a broad wind component 
has been reported before from Ha stationary line analysi s but with 
FWHM reaching only up to 800 kms~^ (Falomo et al. 1987h and 
700 km during the quiescent period reported bv lBlundell et al.l 
(I2OO8I) . although broader widths have been observed before and 
during a flare by Blundell et al. in prep. 

The inner set of narrow lines, moving at speeds ^ 200 km s~ ^ , 
are fairly steady in wavelength, in excellent agreement wi th th e 
presence of a circumbinary ring reported by Blundell et al. ( 20081) . 

The most striking discovery that arises from the Br7 line fit- 
ting is the presence of a pair of widely- separated, hence rapidly 
rotating, narrow components. An example of these narrow (but 
widely separated) pairs are depicted in red and blue colours in the 
upper panel of Fig. |2] The speed with which the radiating material 
spirals in the accretion disc corresponds to half the difference of the 
speed of those lines, under the assumption that the fitted centroids 
correspond to the tangent speed. This reveals material that is spi- 
ralling in the potential well at speeds of about 500 km (see up- 
per panel in Fig.[3])- The accretion disc lines also show a somewhat 
sinusoidal variation, as can be seen in Fig. [3] with a peak-to-peak 
amplitude of ^ 100 kms~^. 

The sixth component, that only appears between (/)orb = 0.2 
and 0.6, seen in Fig.|2]as a small Gaussian blueward of the position 
of the broad wind, may correspond to much more extended line 
emission that is somehow bonded to the system, hence it moves at 
about the systemic velocity but it does not obviously show any of 
the characteristic periodicities. 
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Figure 3. Rotational (upper) and mean (lower) speeds of the rapidly rotat- 
ing components in the Br7 line. Open triangles correspond to the widely 
(500 kms"-*^) separated disc lines, while solid circles correspond to the 
circumbinary ring lines. The rotational speed is calculated from half the 
difference of the redshifts of the red and blue components in Fig [21 Symbol 
sizes correspond to 1 cr error bars. 



4 DISCUSSION 

4.1 The accretion disc and its persistent appearance in the 
near-IR 

Fleeting glimpses of widely separat ed pairs of lines centr ed on Ha 
have been observed very briefly by iFalomo et al] (Il987h . and im- 
mediately prior to a flare by Blundell et al. in prep, implying a 
high rotation speed. It is remarkable that the lines having a rotation 
speed of 500 km seem to be persistent in near-IR observations, 
rather than fleeting or associated with flares as in the optical: we 
have no reason to believe that any of our near-IR observations were 
during, or just before, a flareE Moreover, other near-IR observa- 
tions, both our own in i7-band observations of Brackett-^ 1-74 /xm 
and archival CGS4 observations, are consistent with the presence 
of a pair of lines whose positions correspond to a rotation speed of 
> 500 kms~^. These are depicted in Fig.|4l For a test particle to 
move with a Keplerian velocity of ^500 k ms~^ around a n object 
with a mass of ^16 Mq (as estimated by lBlundell"et aOSo^ it 
would have to be located at approximately 12 Rq away from the 
central object. This gives us an idea of the extent of the accretion 
disc. 

Several authors have utilised hydrogenic lines t o study accret- 
_systems (e.g., iBandvopadhvav et"aD 1 19971 : IShahbaz et al.1 
From obse rv ations of low-mass X-ray binaries. 



Bandvopadhvav et aP (Il999l) have revealed the presence of 
outflows and winds in some accreting mass objects. We emphasise 
the fact that using the Br7 line as a signature of accretion has to 
be undertaken with care since the Br7 emission is telling us about 
outflow via various means, not merely accretion inflow onto a 
compact object. 

The area of each accretion disc line (i.e., the components mov- 
ing at ^ 500 kms~^) seems to be correlated: the difference be- 
tween the areas of these two lines shows a very clear orbital phase 



^ We checked the radio monitoring records of the Ryle telescope (Guy Poo- 
ley) and the RATAN telescope (Sergei Trushkin). 
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Figure 4. Upper panel: Archival CSG4 spectrum of Br7 stationary line 
taken on 1999 July 25 at orbital phase ^orb = 0-03 as a function of helio- 
centric speed in units of kms~^. Lower panel: Brackett-C spectrum taken 
on the last night of our observing run (2006 September 10) at orbital phase 
4>oYh = 0.1 as a function of heliocentric speed. Both observations corre- 
spond to precessional phase ~0.5. The modified Julian day of each obser- 
vation is specified on each plot. 



dependence (see panel la in Fig.|5]). As we can see in panel 2b of 
Fig.[5]the FWHM of each line does not show periodic variations, it 
is around 300 km with a scatter of ^ 70 kms~^. On the other 
hand, panel 3a of the same figure shows that the HEIGHT of each 
line clearly depends on orbital phase. This behaviour is what we 
would expect if the material emitting these components comes from 
the disc, because as the star, together with the spherically symmet- 
ric wind/corona, transits in front of the disc it partially obscures the 
emitting regions. 

The crucial points of the orbital phase dependence are as fol- 
lows: (i) at orbital phase 0.5, the compact object and its accretion 
disc are closer along our line of sight than the companion star is 
— at this orbital phase there is no asymmetric obscuration so there 
is zero difference in the areas of the red and blue lines, (ii) At all 
other orbital phases, the red line area (and indeed height) exceeds 
that of the blue line. From this we draw two conclusions: first there 
is low level obscuration from the stellar wind at all orbital phases 
except when the disc acts as its own windbreak at orbital phase 0.5. 
(We note that most of the time SS 433 's accretion disc is completely 
obscured in the optical.) At all other phases besides 0.5 there is a 
component of the stellar wind between us and the disc lines that at- 
tenuates the blue lines somewhat more than the red lines. Second, 
we conclude that the sinusoidal variation in area difference does 
indeed tell us about the radial extent of the surface of the star, or of 
its obscuring stellar wind. The maximal difference in area occurs at 
orbital phase 1.0 when the star is closer along our line of sight than 
the accretion disc and compact object: this is the phase when the 
greatest attenuating path length though the star (or its dense stellar 
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Figure 5. Panel la: Difference between the areas of the two widely- 
separated (500 kms~^) components, red- minus blue-shifted (the _y-axis 
is in units of 10"-*^^ W m~^). Notice the orbital phase dependence (the 
dashed line is the sinusoidal variation corresponding to that orbital phase). 
Panel lb: Areas of red and blue accretion disc lines as a function of time. 
Panels 2a and 2b are the same as in the upper panel but for the fwhm (in 
kms"-*^) of each accretion disc line. Panels 3a and 3b are the equivalent 
plots for the heights of the Gaussians fitted to the red and blue compo- 
nents, in units of 1.6x 10"^"^ W m~^ /im~^. Error bars correspond to 1 cr 
uncertainties. 



wind) is presented from the accretion disc along the line of sight to 
Earth. The fact that the variation is sinusoidal is consistent with a 
reasonably circularly (spherically) symmetric stellar wind. 



4. 7. 7 Implications for the geometry of the system 

In the previous sub- section we found that the companion star (and 
its wind) eclipses out some areas of the accretion disc. If the orbit 
of the system is essentially circular then the variation of the differ- 
ence of the accretion disc's lines (Fig. [5]) implies that the star has 
a diameter comparable with the radius of the orbit. On the other 
hand, if the orbit is rather elliptical, then the star is smaller and the 
semi-major axis of the orbit must be along our line-of- sight as in- 
ferred from the shape of the light curve, namely, the fact that the 
secondary minima occur s precisely close to mid- way between th e 
primary minima (see eg., iKemp et al.l 19861 : ICroranskii et al.l 1998b . 

Recently, Blundell et al. (200i) found that the total mass of 
SS 433 is approximately 40 Mq. We can roughly estimate the size 
of the system to be about twice the semi-major axis of the orbit a, as 
given by Kepler's third law: — G'Msys(7^orb/27r)^, where Porb 
is the orbital period and Msys is the total mass of the binary. This 
relation implies a size for the whole system of 2a ^ 160 Rq. The 
radius of the companion star should R ^ a/ 2. Thus, the eclips- 
ing of the widely- separated accretion disc lines implies a maxi- 
mum radius for the companion star 7^ ^ 40 Rq. For other re- 
cent mass estimates we obtain similarly large values for the size of 
the star: Hillwig & Gies ( 2008j) report ed Msys ~ 16 M© implying 
R 29 Rq and Ilopcz et aD ( l2006h estimated Msys ~ 55 Mq 
which yields 7? ~ 44 Rq . All of these mass estimates assume the 
case where the orbit is essentially circular. This size is smaller but 
still comparable with the radius of evolved stars undergoing severe 



mass-loss such as 77 Car (^ 80 Rq). This large size is reduced only 
at the expense of accepting an eccentric orbit. The size of the com- 
panion star in the eccentric case would be 6 = Vl — 40 7^© 
where h is the semi-minor axis of the orbit and e is the eccentricity. 

4.2 The disc wind: near-IR view and mass-loss rate 

4.2.7 Precessional- and orbital-phase dependence, or lack 
thereof 

Fig.|2]shows that the br oad wind component , inferred (on the basis 
of optical observations, iBlundell et al .112008b to originate from the 
disc, lacks indications of orbital and nutational periodicities mani- 
fested by the accretion disc. We have not been able to test for pre- 
cessional phase dependence due to our observations spanning only 
0.07 of a precessional phase. A possible explanation of this could 
be that the near-IR wind comes from very close into the centre of 
the disc, where the temperature of the gas could reach ^10^ K. At 
these small radii the material only experiences the influence of the 
compact object and any tidal torque from the companion star could 
be a negligible perturbation. At greater radii as observable in the 
optical, the torques exerted by the companion are much stronger, 
hence precessional periodicities are manifested. 

4.2.2 Mass-loss via line-driven winds 

Close binary systems undergo severe mass loss and mass exchange 
during their lifetime. Like early- type stars, their spectra show a 
number of emission lines due to the large amount of ionised gas 
present in the system. This ionised gas produces thermal radio 
emission and an excess of continuum photons at infra-red wave- 
lengths, which somewhat correlates with line emission, yielding 
similar mass loss estimates. These three spect ral features yield 
mass- loss rate estimates for emission line stars (iKogure & Leund 
I2OO7 



iLeitherej ((19881) has used the luminosity of the optically-thin 
Ha line as a mass-loss tracer. The relation between M and lumi- 
nosity relies on the dependence of both quantities on the density 
of the wind, p(r). Following Leitherer's approach we can derive 
an equivalent expression for Br7 emission, by assuming that the 
emission is optically thin in an isothermal wind. Therefore, the lu- 
minosity of a recombination line, Lu will be given by: 



Ll = Cr^f{T,^)I{vo/Voo,|3) 



R' 



(1) 



where ^;oo is the terminal speed reached by the wind, /i is the mean 
atomic weight and R is the radius of the source (i.e., where the 
wind starts), M is the mass-loss rate, / is an integral defined below 
and (3 i^di parameter that controls the shape of the velocity law of 
the wind. Cn and f(Te^) are constants for each transition given by: 



and, 

/(Teff) 



n{Ai2hvi2 
47rm^ 



2nmek 



3/2 



-3/2gXr./fcTe 



(2) 



(3) 



where ni and 712 are the principal quantum numbers corresponding 
to the upper and lower levels in the atom, respectively. The parame- 
ter bn denotes the departure coefficient and it accounts for non-LTE 
effects (assumed to be ^ 1 for simplicity). A12 is the transition 
probability, k is the Boltzmann constant and Xn is the ionisation 
potential of the level ni. We assume that the wind regions are in 
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Figure 6. Diagram that shows the 2 possible extreme scenarios concerning the configuration of the SS 433 binary system: (left) eccentric orbit, (right) circular 
orbit. The semi-major axis of the orbit is a ~ 80 i?0(see Section U.l.lb . 



radiative equilibrium, i.e., the system's effective temperature Teff 
is very close to the electron temperature, Te. 

The quantity / = I(vo/voo, in equation [T] represents the 
integral over the ratio between the dilution factor and the square 
of the velocity law v{r). In the case where the ve l ocity law cor- 
responds to a /3-law proposed bv ICastor & Lame rs* (1979) andthe 
dilution factor is the standard formula given by Ko gure & Leund 
(l2007h . the integral is given by: 

7=1 r ^-^^ ,d., 

where x = R/r,r being the distance from the centre of the source, 
and u = vo/ Voo ■ The speed vo is the initial velocity of the wind at 
r = R. 

4.2.3 Mass-loss rate via the disc wind 

During our observations, the wind Br7 luminosity attains a maxi- 
mum at 0orb = 0.5, supporting the idea that this emission comes 
from the accretion disc; we have measured this flux density to be 
F(Br7) = (6 it 0.5) x 10"^^ W m'^ The distance to SS 433 
has been accurately determined by iBlundell & Bowlej (l2004h as 
d = 5.5 ± 0.2 kpc. This yields a luminosity for the wind in Br7 of 
14 lb 1 Lq (corrected for = 7.8 magnitudes of optical extinc- 
tion). We expect the wind to start at least at a radius = 20 Rq 
from the compact object, according to the dimensions of the system 
(see Fig. [6]). We computed the integral in equation]?] assuming an 
exponent (3 ^ 0.7 and an initial velocity of the order of the sound 
speed — 15 kms~^ (or 0. 01 Vpo), which is a good approxi- 
mation for hot early-type stars (lLeitherei]|l988h . Because of severe 
dust extinction towards SS 433 there are no measurements of the 
terminal velocity from UV P-Cygni profiles. Therefore we use the 



width of the wind component FWHM ^ 1500 kms~^ as the termi- 
nal velocity of the win d. This is a very m odest value compared with 
some early-type stars (lLeitherejll9"88h . The effective temperature 
has been estimated to be Tgff = 5 x 10^ K by fitting the optical- 
UV spectral flux distribution jGies et al ] |2002l). The mean atomic 
weight of evolved massive stars is /i ~ 2 jpuchs et al ]l2006l) . 

If we apply equation [T]to the Br7 transition of hydrogen (i.e., 
ni — 7 and 712 — 4) we find a useful formula that relates the mass 
loss with the luminosity of the wind: 

Mwind ^ 1 X 10-^ Mo yr-' x /is ^^isoo i?2o^^}o^ (5) 

where /i2 — /i/2, ^;i5oo — ^;oo/1500 kms"-^, i?2o — R/20 Rq 
andLio = L(Br7)/lOL0. 

This value i s the same order o f magnitude as the mass-loss 
rate estimated bv lFuchs et al.l (l2006h who fitted the mid-IR contin- 
uum. These estimates, although discrepant by a factor of 2-3, are 
in resonably good agreement given that both approaches may suffer 
from the assumptions made regarding the geometry of the system 
and the choice of some physical parameters. A crucial unknown is 
the mean atomic w eight of the wind , since the mass-loss ratio scales 
with /i as M oc /i. iLeitherer et al ] ([19951) studied the composition 
of early-type stars and from the H/He ratios they inferred values 
for /i between 1.5 and 2.7. Another important source of uncertainty 
is the choice of velocity law. For /3 < 1 the mass-loss rate scales 
at most by a factor of 2 and it is quite insensitive to the choice of 
vq. On the other hand, for values of p greater than one, M depends 
heavily upon vq. Fortunately, f3 0.7 — 0.8 for winds in 0-stars 
which may match SS 433's hot environment better than B or colder 
stars (which have P > 1). The final important effe ct that must be 
taken into account is the dumpiness of the wind. iLenorzer et al.l 
(I2OO4I) show that the Br7 line may be especially affected by inho- 
mogeneities in the wind. An unclumped wind density can be taken 
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as the real wind density multiplied by a dumpiness factor /. Thus, 
whe n dumpiness is prese nt in the wind M is reduced by a factor 
V7 (iLenorzer et al.lEo04l) . However, the degree of dumpiness is 
not very high closer to the source where the Br7 line is formed, 
thus it should not dramatically affect our estimates. 

Although con tinuum driven w inds might contribute to 
SS433's mass loss (IClarketal.ll2007l) . they are likely to be of less 
importance than line-driven winds in this particular case, because 
of the presence of high velocity outflows (> 1400 kms~^ see 
Section[3l) with high temperatures, diagnostics which rule out con- 
tinuum driven winds as main responsible of mass loss. 



4.3 The circumbinary ring: near-IR view and mass-loss rate 

The idea o f the presence o f a circumbinary disc in S S 433 was con- 
sidere d by iFabrikal (Il993h and recently detected by iBlundell et al.l 
(l2008h by decomposing the Ha line profile in three components 
(during its quiescent state). In the optical this circumbinary ring 
of material orbits the system at a speed of 200 kms~^. It 
has been proposed that the origin of this glowing material cor- 
responds to overflow of gas from the L2 point, assuming that 
the system has filled its Roche lobe (lFabrikalll993l : iBlundell et all 
I2OO8: Filippenko et al. 1988). Our Br7 decomposition is in excel- 
lent agreement with the presence of this excretion disc since we see 
a pair of lines, one permanently blueshifted and the other perma- 
nently redshifted at fairly stable wavelengths, lacking any orbital 
phase dependence and rotating at speeds close to 200 km s~^. 



4.4 Mass-transfer rates in SS 433 and comparisons with the 
Eddington rate 

According to the basic theory of magneto-hydrodynamic winds 
from accretion discs, the mass-loss in the wind is about an or- 
der of magnitude smaller than the mass accreti on rate (inflow at 
the ce ntre of the disc), i.e.. Mace ^10 Mwind (iKonigl & Pudritzl 
Based on our estimate of the mass outflow in the wind 
Mwind = 1 X 10~^ Mq yr~^, this implies an accretion rate of the 
order of Mace ~ 10""^ M0yr~^. On the other hand, the max- 
imum accretion rate allowed by the Eddington criterion, which 
assumes spherical symmetry, is Mace = ^ncrripRxso/ err, where 
Riso = QGM/c^ is the radius of the last stable orbit. This yields 
a maximum accretion rate of Mace ^ 10~^ M© yr~^ for a black- 
hole of 16 Mq. This is greatly exceeded by our estimate of the 
accretion rate in the disc, hence confirming the super-Eddington 
nature of SS 433, which is permitted of course because of the non- 
spherical geometry of the accretion. 



5 CONCLUDING REMARKS 

A succession of near-IR spectroscopic observations, taken nightly 
throughout an entire cycle of SS433's orbit with UIST on UKIRT, 
has revealed: 

(i) The persistent signature of SS433's accretion disc as two 
emission components, blue- and red-shifted, having a rotation 
speed of ^ 500 kms~^. The difference in area and in height be- 
tween each component shows a clear orbital phase dependence, 
very likely the signature of the disc being eclipsed out as the star 
transits in front of it. This contributes only about 15 ± 5 per cent of 
the total Br7 emission line in the near-IR. 



(ii) The presence of the circumbinary di sc (excretion disc) re - 
cently discovered at optical wavelengths bv IBlundell et al.l (l2008l) . 
It contributes about 15 ± 5 per cent of the total emission in Br7. 

(iii) A much faster outflow than has previously been measured 
from the wind. What supports the model that the wind comes from 
the accretion disc is that the luminosity of the broad component 
attains its maximum at (/)orb ^ 0.5, when the compact object is 
closer to us (see Section liJl . The presence of this faster outflow 
{v ^ 1500 kms~^) yields a new upper limit for the mass-loss of 
Mwind = 1 X 10~^ Mq yr~^. This outflow corresponds to a very 
significant 70 ± 5 per cent of the total Br7 line emission. 

These, together with the newly (upwardly) determined masses 
for the components of the SS 433 system, result in an accurate di- 
agnosis of the extent to which SS 433 has super-Eddington flows: 
its accretion rate is 10^ times higher than the Eddington limit. Our 
relatively high spectral resolution at these near-IR wavelengths has 
enabled us to deconstruct the different components, and their phys- 
ical origins, that comprise the Brackett-7 line in this binary sys- 
tem. With this line dominating throughout our series of observa- 
tions by the disc wind, and the accretion disc being only a minority 
(< 15 per cent) contribution, we caution against use of the unre- 
solved Brackett-7 line intensity as an "accretion signature" in X- 
ray binaries or microquasars in any quantitative way. 
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